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excitations of the nucleon. Assuming an effective two-body interaction as 



CO 
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In this paper we develop a theoretical framework which allows us to study 

CN 

in 
o 

\& • a model for low-energy QCD, we derive a relativistic TDHF equation for 

c| ■ a many-body system of quarks. To render the Dirac-sea contribution to 

a: 

the mean field finite, we introduce a symmetry conserving regularization 
■ scheme. In the small amplitude limit we derive an RPA equation. The 

: ^ ■ structure of the ph interaction and modifications due to the regularization 

X 

scheme are discussed. We give a prescription to obtain a nucleon state 
with good angular momentum (J) and isospin (T) quantum numbers on 
mean-field level. To study excitations, we develop a tensor-RPA approach, 
which is an extension of the conventional RPA techniques to systems with a 
nonscalar ground state. This allows us to construct excited states with good 
(J/T) quantum numbers. We discuss a method to reduce the overcomplete 
p/i-space and compute the tensor-RPA interaction matrix elements. Finally 
we extend our scheme to include (| + , |)-states. 
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PACS: 12.39.Ki, 12.40.Yx 

I. INTRODUCTION 

While lattice QCD is on the verge of making definite predictions for certain quantities 
like the static qq potential or the order of the chiral phase transition, hadron spectroscopy 
still remains a difficult task due to finite size effects and uncertainties connected with the 
valence approximation [TJ. Therefore the need for effective models of low energy QCD 
persists, especially in the hadronic sector. Over the years, such approaches have led to 
a deeper insight into various aspects of the relation between the hadronic world and the 
underlying field theory, which is formulated with quark and gluon degrees of freedom. 
We want to focus our attention on models with quark degrees of freedom only, thus in- 
corporating effects originating from instantons or monopole condensation at best in an 
indirect way. Such fermionic models have been quite successful from a phenomenological 
point of view 0. 

For the baryonic spectrum the nonrelativistic quark model, based on the work of Isgur 
and Karl [Q, and related models provide a satisfactory description of the data. Recently 
interest in the nonrelativistic quark model was renewed when it was discovered that a 
spin-flavor interaction allows for the correct ordering of the negative parity states with 
respect to the positive parity ones [|J. 

The mesonic part of the spectrum, at least in the pseudoscalar and, with some restric- 
tion, in the vector channel, was increasingly well understood due to work based on the 
Nambu- Jona-Lasinio (N JL) model |||6j . In that regard the role of the pion as the Gold- 
stone boson of the spontaneously broken chiral symmetry was clarified. The NJL model 
possesses solitonic solutions in the B=l sector which can be viewed as a mean- field de- 
scription of the ground-state baryons. On this basis the phenomenology of ground-state 
baryons could be reproduced quite well [|]||. 
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The aim of the present work is to develop a general framework for a field theoretical 
quark model which describes the nucleon ground state within relativistic Hartree-Fock 
theory and allows to study excited states within an RPA approach. A subsequent paper 
is devoted to the application of our method. 

Wheras in the nonrelativistic quark model the baryon wavefunctions are built from three- 
quark configurations exculsively, we also take into account quark- ant iquark admixtures 
to the wavefunctions. The quark- ant iquark components are included as a coherent super- 
position to the excitations of the valence shell. 

In contrast to the bozonization techniques commonly used in the NJL model, we derive 
a Dirac-RPA scheme in close analogy to the approximation schemes known from nonrel- 
ativistic many-body theory. The essence of such an approach is a particle-hole expansion 
around a mean-field ground state. We emphasise that the methods we develop in part I 
of this paper are to a large extent model independent. Only in part II we make a specific 
choice for the interaction. 

Several versions of a Dirac-RPA approach have been developed in the context of nuclear- 
structure physics In such a relativistic RPA approach it is difficult to achieve a 
restoration of symmetries that are broken by the mean-field configuration. We solve this 
problem on a very general basis by introducing a symmmetry-conserving regularization 
scheme. 

The aspect of symmetry is crucial for our approach also in another regard. Rotational 
symmetry of the interaction in coordinate and isospin space leads to the requirement that 
the physical states carry good angular momentum (J) and isospin (T) quantum numbers. 
In most attempts to construct a self-consistent mean-field state for the nucleon the ro- 
tational symmetries are treated in a rather approximate way. Examples are the Skyrme 
model [11, 12 1 and the NJL soliton model J7]|| which are based on a hedgehog symmetry. 
In such a description the full symmetry group is truncated to an invariance under simul- 
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taneous rotations in coordinate and isospin space, leading to conservation of the grand 
spin G = J + T. 

In contrast to that we construct our states as tensors of the full symmetry group 
SU(2)j ® SU{2)t- On mean-field level this is achieved by a simple projection technique. 
On RPA level the situation is more complicated. The constraints set by the invariance 
properties of the underlying quark model force us to introduce an extended RPA scheme, 
the so called tensor RPA. The tensor-RPA method has been developed in the context of 
nuclear structure physics to study excitations of nuclei with an odd number of protons or 
neutrons [DJ. We generalize this method and adjust it to the Dirac-RPA theory for the 
nucleon. As a result we obtain a powerful scheme which allows us to compute wavefunc- 
tions and excitation energies for states in any (J/T) -channel. 

The tensor RPA is also a quite efficient method to build the correct representations of the 
symmetric group, since it is based on second quantization. In the nonrelativistic quark 
model a complicated classification scheme has to be introduced to construct representa- 
tions of the various symmetries (see, e. g. |14||). Compared to that the construction of 
states via coupled particle-hole operators constitutes an alternative and more straight- 
foreward scheme for building the space of ( J/T)-projected baryon wavef unctions. 
Thus the focus of this theoretical part is on two novel techniques we were forced to intro- 
duce in order to maintain the symmetry structure of a quark model for the nucleon: on 
the one hand a special regularization scheme which guarantees the restoration of broken 
symmetries on RPA level, on the other hand a tensor RPA scheme which allows for the 
construction of states with good angular momentum and isospin quantum numbers. 
The paper is divided into the following sections: In section 2 we derive an equation of mo- 
tion for the one-body density matrix of the system, starting from the Martin-Schwinger 
hierarchy for time-ordered Green's functions. In section 3 we discuss the issue of con- 
sistent implementation of a regularization scheme. We focus on the Dirac Hartree-Fock 
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problem for the nucleon in section 4. Section 5 is devoted to the discussion of excited 
states belonging to multiplets of the group SU(2)j x SU{2)t- The paper closes with a 
discussion and a summary of the results. 



II. EQUATIONS OF MOTION FOR A RELATIVISTIC SYSTEM OF 

FERMIONS 

In the framework of field theoretical quark models the common way to derive many- 
body equations of motion is via bosonization techniques [HJ. There one considers an 
Euclidean path integral and derives the equivalent of the Hartree-Fock approximation by 
restricting the functional integration to the stationary phase configuration. The equivalent 
of RPA can be derived by allowing for a time dependence of the mean field. 
We follow a different path, because we want to stay as close as possible to the usual 
formulation of many-body physics in terms of multi-particle multi-hole operators acting 
on a mean-field ground state. 

Since we do not perform a Wick rotation to Euclidean space-time we avoid ambiguities 
with the treatment of vector potentials [jTJ. We start from a model state given by the Dirac 
sea with three valence quarks on top. Dynamic correlations are treated on the level of 
small-amplitude oscillations of the one-body density matrix. In that spirit we derive many- 
body equations of motion using a relativistic version of the Martin-Schwinger hierarchy 
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The class of models we consider is defined by a Lagrangian of the form 

L = J dW(l)(ij^ - m)*(l) --J dld2dl'd2'(12\u\l'2')^(l)^(2)^(2')^(l'), (1) 

where 1, 2, 1', 2' denote the quantum numbers necessary to label the single particle states, 
including the time label. It is assumed that the quarks interact via an effective two-body 
potential u. We define time ordered Green's functions by 



G(l . . . n, 1' . . . n') = (T[*(l) • • • *(n)*(nO • • • (2) 
With that we obtain an equation of motion for the two-point function 

D(1)G(1, 1') = 5(1 d2dl"d2"<12|tt|l"2")G(l"2", 1'2+) (3) 

and the adjoint 

L>*(l')G(l, 1') = 5(1 -!')-«/ rf2dl"d2"(l"2'>|l / 2)G(12-, l"2"). (4) 
The Dirac operator £) is given by 

D = ry M <9 M - m (5) 

and the notation _D(1)G(1, 1') denotes the action of D on the first argument of the Green's 
function. The superscripts +(— ) indicate that the argument has to be evaluated at an 
infinitesimally later (earlier) time compared to the integration variable. 



Following closely the derivation of time dependent Hartree-Fock theory given in |L8 
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approximate the 4-point function by an antisymmetrized product of 2-point functions 

G(12, 1'2') = G(l, l')G(2, 2') - G(l, 2')G(2, 1'). (6) 
Furthermore we introduce a one-body density matrix by 

p(l,l')=^(l,l')7o(l'), (7) 

where 70 is the Dirac matrix times a unit matrix with respect to space-time, color and 
flavor indices. 

If u parametrizes a one-boson exchange interaction with a local vertex we encounter a 
divergence in the equations of motion, because in this case one of the density operators 
appearing in the integral is evaluated at the same points in space-time. Thus, the theory 
requires regularization. In section 3 we will discuss a regularization scheme which strictly 



preserves the symmetries of the equation of motion. For the time being we perform all 
algebraic manipulations as they can be applied to finite quantities. 
We introduce new time variables by 

T = h-t v , t=(t 1 +t r )/2. (8) 

Multiplying Eq. (||) and @ by 70 and taking the difference leads to an equation which 
only contains a time- derivative with respect to t 

[id t - (a ■ p + (3m)(l) + (a ■ p + (3m)(l')} p(l, 1') 
= - J d2dl"d2"(l2\v\l"2")p(l",l')p(2",2 + ) 

+ J d2dl"d2"(12\v\l"2")p(l,l")p(2-,2"). (9) 

In this expression 

v = 70 <S> 70 u (10) 

plays the role of a two-body potential as defined in a hamiltonian formulation. The 
shorthand v indicates antisymmetrization with respect to particles 1 and 2. The tensor 
product of the two Dirac matrices defines a two-body operator. In the above derivation 
we have assumed that v is hermitian. 

If the interaction is local in time, we obtain a closed equation for the evolution in t 

idtp(l,l',t) = [H MF ,p](l,l',t), (11) 
where the limit r = t% — ty — > 0~ was taken. The mean-field Hamiltonian is defined by 

H MF (l,l',t) = [a •p + /3m](l,l')+ / d2d2'(12\v\l'2')p(2', 2, t). (12) 
We thus recover the same structure as in the TDHF equation known from nonrelativistic 



physics [19]. Although not written in a manifestly covariant manner, Eq. (|TT|) is Poincare- 



invariant if the same is true for the two-body interaction v. This can be traced back to 



the approximation introduced in Eq. @, which respects the covariance properties of the 
4-point function. 

In most casesQ Eq. (|TT| ) is of little use for practical purposes unless we specify a regu- 
larization scheme. The two-point function will be infinite due to contributions from the 
Dirac sea, so we need to filter the spectral density of negative energy states to obtain 
finite results. 

In the next chapter we will discuss a method that allows us to regularize the TDHF- 
equation in a consistent way. This method will lead to the derivation of a self-consistent 
Dirac Hartree-Fock + RPA approach. 

III. REGULARIZATION OF THE TDHF EQUATION 

A covariant theory as the one developed in the previous chapter confronts us with 
single particle spectra which are unbound from below. In that case the standard RPA 
approach for the excitations of the system proves to be inappropriate. Even when one 
considers a mean field that is made up by positive energy states only, one has to take 
special care to develop a consistent RPA scheme, because the response function picks up 
contributions from negative energy states |§. 

We use the term 'consistent' in the following way: Given a Hartree-Fock ground state 
a method to compute the low-lying excitations of system is consistent if it respects the 
symmetries of the underlying theory. If a symmetry is broken in the ground state, a 
mode that corresponds to a symmetry transformation (e. g. a translation in space) has 
to appear at zero excitation energy. Then the spurious modes separate from the intrinsic 
excitations. A scheme that does not preserve the symmetries will lead to unphysical ad- 
mixtures of spurious modes to the excited states. 



1 There are, however, a few exceptional cases where no regularization is needed (see, e. g. pOfl). 
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In the nonrelativistic case the consistency problem has been discussed in much detail 



2l| , |22[| . The most important result is that RPA modes corresponding to symmetry trans- 
formations appear at zero excitation energy. However, a straightforward extension of the 
nonrelativistic RPA scheme does not suffice to describe the linear response of the Dirac 
sea. Such a scheme does not allow for the decoupling of spurious modes in the relativistic 
case. 

If the mean field is exclusively made up from positive energy states the simplest way 
to achieve consistency is to shift the negative energy poles of the Feynman propagator 
into the lower half plane [[|. This corresponds to a drastic change of the structure of 
Dirac-hole states from backward to forward motion in time. Effectively this simulates an 
unoccupied Dirac sea. As a consequence the RPA equations include contributions from 
configurations where a particle is scattered into a Dirac-hole state. 

We follow a different path for the regularization of our time dependent mean-field theory. 
On the one hand we include explicitly the effects of an occupied Dirac sea. On the other 
hand we avoid the appearance of unphysical states in the ph basis. On the RPA level our 
regularization scheme amounts to using a modified (with respect to the original theory) 
interaction. The eigenstates are defined as a superposition of the physical (positive en- 
ergy) ph states and their time-reversed analogues. 

The basic idea of our regularization scheme is to replace the divergent mean-field Hamil- 
tonian, which is a functional of the one-body density, by an effective Hamiltonian of the 
form 

H' MF (p) = H MF {R[H' MF )pR[H' MF \). (13) 
Thus only contributions to the mean field from the regularized density 

p' = R{H' MF ]pR[H' MF ] (14) 
are taken into account. R[x] is a cut-off function with the property 



R[x] -»• 0, (z -> -oo) (15) 

for a real argument x. The falloff for large negative arguments should be fast enough to 
guarantee that the effective particle number of the ground state J2iR 2 [^i] remains finite. 
The sum is to be taken over the occupied single-particle states with energies e^. We 
assume that the cut-off function can be expanded into a power series, which provides the 
definition of R[x] for an operator-valued argument. Written down in a specific model space 
Eq. fljjD is an implicit relation between the matrix elements of the effective Hamiltonian 
H' MF , which appear on both sides of the equation. 
The modified TDHF equation reads 

id t p=[H' MF (p),p\. (16) 

The new description allows to solve for a set of states that simultaneously diagonalize p 
and H' MF , corresponding to the stationary limit of Eq. When we write down H' MF 

in a basis of such eigenstates, contributions from states arbitrarily deep in the Dirac sea 
are excluded by the cut-off function R. However, the density operator itself will still 
diverge in a basis of eigenstates of the position operator, corresponding to a completely 
filled Dirac sea. 

The regularization scheme we propose preserves the symmetries of the underlying inter- 
action v. To show this, consider a unitary transformation 

p(9) = lPp(0)U, (17) 

with 

U = exp(i9 a S a ). (18) 

S a are the generators of any group of transformations that leave the interaction invariant. 
Let us assume that p(0) is a solution of Eq. flTE|). If we can show that H' MF transforms as 
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H' MF {9) = rfH' MF (0)U, (19) 

we have clearly proven the invariance property of the regularized equation of motion. In 
the above expression the abbrevation 

H' MF {9) = H' MF {p{6)) (20) 

is used. 

From the defining equation for H' MF we have 

H' MF {9) = H MF {R[H' MF {e))p{9)R[H' MF {9)]). (21) 

The proof of the transformation law Eq. ( |19"D is based on the fact that Eq. (|2T| ) is to be 
considered as an implicit definition of H' MF (6). We assume that we are given a solution 
H' MF (0) of this equation for 9 = 0. To prove the validity of the transformation law, we 
insert the right hand side of Eq. ( |19|) for H' MF {6) 



U^H' MF (0)U = H MF (R[U^H' MF (0)U] p{0)UR{tf H' MF {0)U]). (22) 

The cut-off function transforms as 

R[U^H' MF (0)U] = U^R[H' MF (0)}U, (23) 

as one finds from a power-series expansion for R. The transformation law for the original 
mean-field Hamiltonian is [ETC 



H MF (rfp(p)U) = U^H MF (p(0))U. (24) 

With that Eq. flU) becomes 

tfH' MF (Q)U = U^H MF (R[H' MF (0)}p(0)R{H' MF (0)})U, (25) 

which is the defining equation for the regularized mean-field Hamiltonian at 9 = 0. Thus 
the expression given in Eq. (|T9"| ) is a solution of the defining equation whenever the same 
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is true for H' MF (0). Note that a proof of the transformation law can only be given in such 
an indirect manner, since Eq. (pl|) can in general not be solved for H' MF (8). 
Our regularization scheme relies on the fact that the cut-off is provided by a description 
which is self-contained, i. e. with no reference to an external set of basis states. Had we 
done an expansion with respect to any complete set of states, e.g. momentum eigenstates, 
we would have violated the symmetry requirements if we had cut off the high momentum 
tails of the density matrix. 

In fact we have introduced a whole set of regularization schemes, because so far the form 
of the cut-off function R has not been specified. Clearly the many-body-theory defined 
above is not unique and in general we expect the results to depend on the cut-off function. 
Any parametrization which includes the effects of the most important states close to the 
fermi surface should, however, lead to similar results for physical observables. This can 
only be checked by applying the relativistic TDHF theory to a specific model. A similar 
situation is encountered in the NJL-soliton model, where different schemes are used to 
regularize the diverging contributions from the Dirac sea. 

We now proceed to develop an RPA equation by considering the small amplitude limit of 
Eq. (|16|) . For this purpose we write the density operator as 

P = Po + Sp. (26) 

Po is time independent and fulfills 

[H' MF (p ),Po] = 0, (27) 

5p is a small perturbation. Inserting this ansatz into the TDHF equation and keeping 
terms linear in Sp we obtain 

l9t5p = [ ^T UP ' Po] + VTurM' 5p] (28) 
In the unregularized RPA we would have obtained an expression where H' MF is simply 
replaced by Hmf- The second term of Eq. ( P8|) contains the information about single- 
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particle energies, since we can make use of the fact that H' MF and po are simultaneously 
diagonal. The derivative in the first term defines the ph interaction which enters the RPA 
equation 

(ijKffli'f) = ^f^k (29) 

A similar definition of an effective interaction enters the Landau-Migdal theory of Fermi 
liquids [|3[ . 



In the standard derivation of the RPA equation one uses the fact that all particle-particle 
and hole- hole matrix elements of bp vanish |2l|] . In case that the mean- field ground state 



can be written as the vacuum of a set of single-particle annihilation operators the same 
statement will hold for the Dirac RPA. 

The nucleon, which is a system built from a filled Dirac sea with three valence quarks on 
top, requires some special care. If we model the nucleon by single-particle wavefunctions 
taken from a mean-field calculation, we have to introduce explicit correlations to construct 
a (J=l/2,T=l/2) ground state. In a strict sense this does not allow us anymore to define 
the terms 'particle' and 'hole'. However, the old concepts will remain unaltered for all 
states outside of the valence shell. Since such ground-state correlations are restricted to 
the valence shell, we do not expect to make a big mistake if for the derivation of the RPA 
equation we ignore them. 

A similar situation is encountered in nuclei away from magic numbers. Already by statis- 
tical reasons such systems are characterized by ground-state occupation probabilites less 
than unity in the valence shells. Superimposed are dynamical correlations due to residual 
interactions. In the theory of Fermi-systems the various contributions to such ground- 



state correlations are being subsumed into the quasiparticle concept [23J. This picture 
takes into account that the mean-field model states, e. g. from Hartree-Fock theory, are 
distributed over a range of many-body eigenstates. 

We now focus our attention on the effective ph interaction which originates from the reg- 
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ularization of our theory. In order to compute the interaction matrix elements we have 
to examine the derivative of the regularized mean-field Hamiltonian with respect to the 
density. The derivative has to be evaluated at po; the stationary solution of the TDHF 
equation. In the following all single particle labels refer to states of the self-consistent 
basis, i. e. the set of states which simultaneously diagonalize the regularized mean-field 
Hamiltonian and the density matrix. As shown in Appendix A, the implicit equation for 
the effective interaction is obtained as 







-m 



dpkk 
dH MF 



MF 



' IF -R[e k ]R[e k i 



+ 



PO 

d{H' h 



'mf)w 



po d PM> 



U[Cl ' ] " /i>h] (9(k F - l)R[ £l ] + 6{kp - l')R[e v \) . (30) 



po 1 



dp'iv 

As an example which reveals the structure of the above equation let us consider a cut-off 
function of the form 

R l [e i ]=e-*- KH+K \ (31) 

where / has dimension [MeK -1 ]. Clearly Ri has a power series representation and 

limine*] =0(e; + A). (32) 

In the following let us restrict ourselves to very large Z, where the step function is repro- 
duced with high accuracy. We want to emphasize that the above choice for the cut-off 
function is made in order to provide an interpretation of Eq. ( |3T)D in terms of ph states 
and to compare our approach to the previously mentioned Dirac-RPA method of Dawson 
and Furnstahl |J. Their approach allows for particle states in the Dirac sea. At the 
moment we do not touch the question if the special choice for the cut-off function in 
Eq. ( [31]) is suited for practical purposes, because our aim is to illuminate the structure 
of the equations rather than going into details of a numerical solution. In fact, as has 



been shown in Ref . |25[ , a sharp cut-off might lead to problems when attempting to solve 
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self-consistent mean-field equations by numerical methods. 
We write the matrix elements of the effective interaction as 

rr _ d(H' MF )kk> /qq\ 

Vkk',11' = a • K 66 

Opiv 

When two-body states a = (kk'), (5 = (//') are introduced the defining equation for the 
effective interaction can be written as a matrix equation of the form 

Vcf, = M aix V^ + N a ^ (34) 

where the matrices M and N can be read off from Eq. (|30|). We introduce model spaces 
8 and T> of two-body states defined by 

8: ph states (kk 1 ) with ^ f or [k <-> k') 

— A < 6fc < €f 

V: states (kk') with ^ ^ or (k «-> fc'). 

e*' < -A 

denotes the Fermi energy. Clearly the spaces 8 and T> are orthogonal. 
Let P and Q be the projectors onto £ and P. For an arbitrary operator A we introduce 
the notation App = PAP, Aqq = QAQ, similarly for other combinations of P and Q. 
Note that P and Q are operator- valued functionals of the self-consistent density p . 
The structure of M in Eq. is such that the sum over \x is restricted to states in T>. 



This fact is due to the special choice for the cut-off function Eq. ([H]). With that we 
obtain the two relations 

V PP = MpqVqp + Np P 

V QP = M qq V Q p + N QP . (35) 

When we solve the second equation for Vqp and insert the solution into the first equation, 
we find 
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Vpp = M PQ (1 - Mqq)- 1 N qp + N PP . (36) 

Similar structures involving projection operators are encountered in the theory of valence 
forces in nuclei [26|. There one considers the modification of the interaction due to the 
truncation of the Hilbert space to a smaller model space. 

For the TDHF equation we exclusively need ph matrix elements of V, i. e. matrix elements 
in S. The model space T> contains the de-excitation states of Dawson and Furnstahl ||, 
which can be viewed as negative energy ph states with respect to a 'Fermi surface' at cut-off 
energy. In contrast to their prescription the de-excitation states only enter as intermediate 
states in our RPA scheme and do not contribute to the collective ph amplitudes describing 
the excited states of the system. 

When the RPA eigenvalue equation is solved with the effective ph interaction defined by 
Eq. fl36|), the spurious modes decouple from the excitation spectrum. 



IV. MEAN-FIELD THEORY OF THE NUCLEON 

For the ground state of the nucleon we have to solve the stationary limit of the reg- 
ularized TDHF equation, which amounts to finding the eigenvalues of the regularized 
mean-field Hamiltonian. The fact that Hmf is replaced by H' MF formally leads to a 
highly nonlinear problem. The eigenvalues are, however, determined iteratively. 
In order to exemplify the procedure without the quite complicated valence-shell contribu- 
tions to the mean field, let us consider the analogous problem for the vacuum. Initially 
we choose a set of states which should not be too far away from the self-consistent set in 
order to guarantee convergence. We then compute the regularized mean-field Hamiltonian 
assuming that the states we have chosen are eigenstates. The Hartree-Fock potential of 
the vacuum acquires the form 

Vmf = ^(pi\v\qi)ala q R 2 [ei}8(-ei), (37) 

pqi 
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where in the argument of the regularization function H' MF has been replaced by an ap- 
proximation for its eigenvalues. After that we solve the eigenvalue problem for the new 
H' MF and repeat the procedure until convergence is achieved. By construction a configu- 
ration obtained in such a way represents a solution of Eq. Q27p. 

A well known model for vacuum structure and spontaneous chiral symmetry breaking is 
the Nambu-Jona-Lasinio(NJL) model (for a review, see e. g. |27|). Due to the simple 
momentum structure of the NJL model it is easy to derive a nonlinear equation which de- 
termines the chiral condensate (\E'\E'). With the assumption of a translationally invariant 
vacuum state this so called gap equation can also be derived within our Dirac-Hartree- 
Fock scheme applied to the NJL model. One finds 

oc / dEVE 2 -m 2 mR 2 [E], (38) 

where m is the constituent quark mass which includes contributions from the scalar part of 
the mean- field potential. Eq. ( p8[ ) is identical to the gap equation with an 0(3)-invariant 
cut-off discussed in [^7[] when a step- function parametrization is chosen for R 2 [E]. 
From this discussion it should become clear that the Dirac-Hartree-Fock method set forth 
in the previous chapter contains previously developed methods for the determination of 
the mean-field structure of fermionic systems. The special form of the regularization 
scheme based on an implicit definition of the mean-field Hamiltonian mainly has an effect 
on the RPA sector of the theory where it leads to a modification of the ph interaction. 
Compared to the vacuum sector additional complications are introduced when one in- 
cludes valence quarks and studies a nucleonic system. As mentioned in the introduction, 
our aim is to construct states that carry good J/T quantum numbers. However, the J/T 
symmetries are not automatically self-consistent symmetries of the nucleon mean-field 
Hamiltonian. This means that the mean-field of a system described by a many-body 
state with good J/T quantum numbers is no longer invariant under rotations in coordi- 
nate and isospin space. Subsequently it will be only be possible to fulfill the symmetry 
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requirements by acting with a projection operator on the eigenstates of this symmetry- 
broken Hamiltonian. 

A possible approach would be to define the nucleon as a state with grand spin G = 0, 
where the grand spin is defined as 

G = J + T. (39) 

This so called 'hedgehog' ansatz is used in the chiral soliton models, because the grand 
spin symmetry is a self-consistent symmetry of the mean-field Hamiltonian. After solving 
the mean-field equation states with good J/T quantum numbers are projected out, usually 
by a cranking procedure. In the language of nonrelativistic nuclear physics this approach 
would be called a 'variation before projection' method |2l] . 



A superior but technically more involved approach is provided by the 'variation after 
projection' techniques. There the mean-field approximation is applied to a projection 
of the full hamiltonian onto a subspace of J/T eigenstates (see f28| , |29"| and references 
therein). The projected ground-state wave function can in general only be written as a 
superposition of several slater determinants. 

In order to keep the theoretical discussion more transparent and numerical investigations 
feasible we restrict ourselves to a somewhat simplified 'variation after projection' technique 
in the following. We assume that the valence shell is a k = 1 or k = —1 state, thus carrying 
an angular momentum of J = 1/2. Although the notation we use suggests a valence level 
of positive energy, the positivity of the valence energy is not mandatory for our discussion. 
The vacuum |0) which we refer to as 'Dirac sea' is defined by the single-particle states 
which are lower in energy than the valence level. 

A wavefunction with the quantum numbers of the nucleon ((J P ,T) = (| + , |)) is given by 



l r "> r *> = 4f E T &LrAlAlAlM- (40) 

V 10 r 2 i,r z2 ,r z3 
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The quark states are labelled by their color quantum number (superscripts) and angular- 
mo mentum/isospin quantum number (subscripts). The (J,T) quantum numbers of the 
nucleon are denoted with = (1/2,1/2), the projection quantum numbers with T z . 
In the following we take the convention that doublets of (J/T) quantum numbers are 
subsumed in one single symbol 

V^iJ^Ti). (41) 
Projection quantum numbers are abbreviated in a similar way 

r zt = ((J z ) t ,(T z ) t ). (42) 
Unless stated otherwise, we define for any function of the (J/T) doublets 

T 2 , . . . , r n ) EE /( J X , J 2 , . . . , J n )/m, T 2 , . . . , T n ). (43) 

State vectors and tensor operators which are labelled by T-symbols are an exception to 
this convention. To express the nucleon wavefunction in terms of creation operators we 
have introduced coupling coefficients ' r -p of the group SU(2)j <g> SU(2)t @. The 
T-coefficients do not split up according to Eq. (|4"3"D. 

For the further discussion it is appropriate to introduce a creation operator for a 3-quark 
state with nucleonic quantum numbers 



By definition this operator fulfills 



1 ^ r 2l ,r z2 ,r z3 



Ai(r z )\o) = \r N ,r z ) 

A N (T z )\0) = 0. (45) 

The diagonal matrix elements of the density operator defined with respect to a nucleon 
state with projection quantum numbers T z then become 
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p V*\m) = (0\A N (T g )*\m)*(m)A]f(r t )\0) 

= (o|^( m )^( m )| ) + (o|[[A iV (r 2 ) ) *t( m )^( TO )] ) A] v (r z )]|o}. (46) 

From this expression it becomes obvious that the density operator contains contributions 
from the Dirac sea and from the valence quark configuration, respectively. It is easy to 
verify Eq. fl4"5|) by writing out the double commutator and making use of Eq. (^Sj). The 
different contributions are evaluated as 



p(r.)( m ) = g( r *\ m )/18 + 6(-m), (47) 

where the symbolic 9 function is defined as 

{1 : m e Dirac sea 
(48) 
: otherwise 

The valence-occupation factors g are color independent and given by 

J2r 1 r n(Tr* r r J 2 : m £ valence shell 
g (^\m) = J ^ r ^ K ^r zl r z2 J (4g) 

: otherwise 

From these relations it is apparent that the valence part of the full wavefunctions leads 
to a dependence of p^ z \m) on the J/T projection quantum numbers of m. In a basis 
of eigenstates of H' MF the mean-field potential for a nucleon with projection quantum 
numbers T z is defined as 

Vmf = E(P^I^)4 a ^ 2 hlP (r2) W- (50) 

pqi 

Due to the noninvariant structure of the density matrix the mean-field potential is no 
longer invariant under rotations in coordinate and isospin space. Such complications 
are well known from nuclear structure physics and have been treated on mean-field level 



allowing for angular momentum and isospin mixing in the wavefunctions [31 



In the present context the approximation scheme we use to describe the excited states 
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requires a single particle basis of J/T eigenstates. In order to enforce that the mean-field 
potential is invariant with respect to rotations in coordinate and isospin space we use 
occupation probabilities averaged over projection quantum numbers 

This approximation neglects the nonscalar parts of the density matrix. As a result, the 
averaged density becomes independent of the projection quantum numbers and is given 
by 

p{m) = - A (52) 

for any state m of the valence shell. 

Comparing this prescription to the usual hedgehog ansatz we find that p contains contri- 
butions from configurations with G / 0. Due to the unconstraint summation over angular 
momentum and isospin projection quantum numbers G z can take the values 0, ±1. 
Another difference with respect to the hedgehog ansatz is related to the Lorentz-structure 
of the mean field. Due to the averaging procedure only the Lorentz-scalar and the time-like 
component of the Lorentz-vector part contributes. Pionic contributions to the stationary 
mean field which are a major building block of the hedgehog soliton are thus excluded. 
Note, however, that before the averaging the expectation value (TN^z^lb^^N^z) 
can take a nonzero value. For the same reason the pion field will contribute in the more 



elaborate projection scheme of Refs. [28,29|. With such a description the tensor- RPA 



matrixelements to be discussed in the following section will exhibit a quite complicated 
structure. Thus in the present work we restrict ourselves to the averaging prescription of 



Eq. ([51]). 

Within this simplified scheme the pionic contributions are recovered on the level of ground- 
state correlations. Such correlations are inherent in an RPA description, where the in- 
formation on improvements of the ground state beyond the static mean-field level enters 
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via the ^-amplitudes. In [^2| a simple RPA scheme for the nucleon was proposed where 
the tensor structure of the excited states is ignored. In such a scheme the single particle 
occupation probabilities for the RPA ground state are given by 

n h = 1 - (2T h + l)- 1 \ £ (2T + l)\Y; h (T,T z np(h)-p(p)y 

n p = (2T h + iy l1 - £ {2Y + l)\Y p \{T,T z )\\p{h)-p{p))\ (53) 
z r,r*,v 

where rih and n p are the occupation probabilities for hole and particle states, respectively. 
The backward amplitudes Y are defined in the usual way (see, e. g. ||21|| ). The depletion 
of hole states and occupation of particle states is due to 2p-2/i-admixtures to the ground 
state. These admixtures contain contributions from ph excitations of the vacuum, where 



the pionic modes are expected to be most important. In |J2| it has been shown that the 
vacuum contributions are in fact sizeable, leading to a depletion of the Dirac sea of about 
10%. 

With such an improved description of the ground state not only the coupling to pionic 
modes is recovered, but it is also possible to study the spin structure of the nucleon. In 
the nucleon wavefunction of Eq. fl40|) the spin is carried by the valence quarks exclusively. 
When ground-state correlations are taken into account, a part of the angular momentum 
is transferred to ph admixtures. 



The ground state density of Eq. (f£6|) or Eq. fl5T|) does not allow an unambiguous assign- 
ment of particles and holes, since the valence shell is only partially filled and the vacancy 
can be populated in excitations. The structure of the ph operators to be included in the 
RPA basis will be discussed in the next section. 

Applying Eq. ([51]) to a valence quark configuration coupled to Ta = (3/2,3/2) the same 
averaged occupation probabilities as in Eq. (|52|) are found. Thus, the (l/2 + , 1/2) and the 
(3/2 + , 3/2) valence quark configuration are indistinguishable on the level of the averaged 
density matrix. However, as discussed in the next section, this degeneracy will be lifted 
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from configuration mixing from residual interactions. When we take the expectation value 
of the full Hamiltonian with an explicit representation of the ground state like Eq. ( f40"|) we 
in general find different values for nucleon and A. This means that the expectation value 
of the residual interaction is nonzero for both states due to ground-state correlations. 
For the derivation of the RPA equation in the previous chapter we have assumed that 
these correlations are neglible, i. e. we have derived the equation in the same manner as 
for a closed-shell system. A measure for the importance of correlations is the difference 
|(rjv|-ff|r/v) — (rA|-£f |Ta)|. If this difference is not small compared to the lowest single 
particle excitation energies a RPA scheme assuming an uncorrelated ground state must 
be taken with precaution. 

The important point of our description is the fact that we can extract a single-particle 
basis with good J/T quantum numbers from the mean- field Hamiltonian defined with the 
averaged density. This will be a crucial prerequisite for the tensor-RPA method discussed 
in the next section. 



V. TENSOR RPA 

The tensor-RPA approach is a special case of the tensor equations of motion method 
developed more than 20 years ago by Rowe et al. which, however, never seems to have 
been used in realistic calculations. It allows to compute excited states of a system with 
a nonscalar (with respect to a S'?7(2)-symmetry) ground state. Usually RPA calculations 
are restricted to closed-shell systems with a spin/isospin zero ground state. In this case 
the excited states carry the quantum numbers of the ph configuration admixed to the 
ground-state wavefunction. 

For a nonscalar system the situation is more complicated. We assume that the Hamil- 
tonian which governs the dynamics is invariant under spatial rotations (SO (3)) and 
isorotations (577(2)). Since SO (3) and SU(2) are isomorphic, the invariance group is 
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SU(2)j ® SU(2)t- The eigenstates of the Hamiltonian should transform as tensors with 
respect to this symmetry. An excited state of given J/T quantum numbers can in general 
be defined by coupling ph pairs of different tensor ranks on the ground state, according 
to the rules for outer products of SU (2) representations. 

For the nucleon we have to set up such an extended RPA formalism in order to com- 
pute excitations. For the time being we want to stay in the SU(2) sector of the flavor 
group. This seems to be a reasonable approximation, because hints for strangeness in 
the excitation spectrum of the nucleon are restricted to the iV(1535) resonance, which 
decays to a sizeable amount into a nucleon and an rj meson. A SU(3) tensor RPA could 
be formulated as a straightforward extension of the method discussed here, but in order 
to apply it to baryon phenomenology one would have to think of a mechanism which 
splits the multiplets into several SU(2) representations and which can be traced back to 
a symmetry breaking mass matrix ||J]. 

The formal derivation of the tensor-RPA equation for the nucleon is independent of a 
specific representation of the ground state. The general formalism is therefore applicable 
also in the context of a more sophisticated description than the one used in the previous 
chapter. Details of the method as the number of linearly independent basis states or 
the structure of the interaction matrix elements, however, depend on the ground-state 
structure. 

A. The Tensor-RPA Equation 

We assume that an excited state of the nucleon is defined by the action of an excitation 
operator on the ground state 

l*r A » = (Ql x \T N ))) r \ (54) 
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rAtTjv) denotes the SU{2)j ® SU{2)t tensor rank of the excited state (ground state) and 
x includes all additional quantum numbers. Since the adjoint of a tensor state is found 



to transform as a tensor only when supplied with appropriate phase factors ||13|| , we use 
the definition 

|r,r,» = |r,r 2 > 

«r,r,| = (-i) r *(r,-iU (55) 

where V z is the projection quantum number. A similar definition is used for operators 

0\T,T z ) = 0\T,T z ) 

o(f,r,) = (-i) r+r *o(r,-r z ). (56) 



The notation for the tensor coupling we use in Eq. (|5^), applied to arbitrary tensors R 
and S, stands for 



(r^ x s r >f r = £ (rarar^irr,)^^, (57) 



r 

r z ir z 2 

with Clebsch-Gordan coefficients that split up according to Eq. ([43]). The excitation 
operator Q^ A contains a sum over several tensor ranks Tj 

Qk = E°iiV (58) 

i 

according to angular-momentum/isospin selection rules. We demand 

Qa\T n )) = (59) 

to ensure that the excited states are orthogonal to the ground state. Furthermore we 
assume that the excitation operators can be written as a superposition of certain operators 
yet to be specified 

oi ri =E%(^( r *)- ( 6 °) 

p 
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Based on these assumptions a general equation of motion for the excited states of the 



nucleon can be derived, which in structure resembles the results of Ref. 13 



MZ{x) =lu a NZ(x), 



(61) 



where are the excitation energies. The hermitian matrices M and N are defined as 



Mai M = Y,n-^ rN - rA ~ r - r '^N\\K(r t ),H, v l(r 3 )] r \\r N ) 

r 

•w^r^r^rTA) 

N ai , Pj = i:r(-i) r -- r -- r - r ^r JV ||[r 7ci (r,),r 7 t(r,)] r ||r JV ) 

r 

■^(r^rvr^rrA). 



(62) 



In this expression H is the full Hamiltonian of the system and W are the SU(2) j®SU(2)t 
Wigner 6j-symbols |33] . We have introduced a frequently needed symbol for the dimension 
of multiplets by 



f = (2r + 1)2. 



(63) 



A coupled commutator of two operators P and Q is defined as 



P v \Q l 



E (rirar^irr,)^,^, 



r 2 ir z 



A symmetrized double commutator will be denoted by 



(64) 



2[A,B,C] = [A, [B,C]] + [[A,B),C] 



(65) 



When one of the three operators is a scalar, as in Eq. (|62|), we use the notation of Eq. 
(|B~4l) for the coupling of the remaining two operators. 

Depending on the choice for the operators r\ one can derive several approximation schemes 
from this equation, e. g. Tamm Dancoff, RPA or quasiparticle RPA. The states are nor- 
malized according to 
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Z\y)NZ(x) = S x . 



xyi 



(66) 



with a metric tensor N which will in general be different from the identity matrix. For 
deriving an RPA equation, we make the specific choice 



O j x r k = £ [X( m i)k{x)A ] mi (T k ) - Y {mi)k (x)A mi (T k ) 



(67) 



An operator A' mi (r k ) creates a ph pair (mi) with J/T quantum numbers T k and is defined 

as 



A-mi^Xk^zk) — X! (^mX'i^zm^zi\^k^zk)0'r m r zrn a rir zi - 



(6* 



We adopt the convention that particle states are denoted by m,n and hole states by 
A state above the valence shell is always a particle, a state in the Dirac sea always a hole. 
Within the valence shell there can be both types of states. This opens up the possibility of 
having ph states where both particle and hole are in the valence shell. Such states would 
induce a recoupling of the ground state and can be used to express the wavefunction of the 
A(1232), which is degenerate with the nucleon on mean-field level. However, introducing 
such recoupling states into the set of excitation operators requires special care. This issue 
will be discussed later. 

With Eq. (p7|) the tensor equation of motion takes the form 





-M( 2 ) 


\ 




I 












v -M< 3 ) 


M< 4 ) 


/ 




V 



AT(l) _AT(2) 
_AT(3) 



(69) 



With the definition 

Cl,i = (-1)^-^-^ WiTkT^^NjTA) 
the submatrices are given by 

C l k \i(T N \\ [A mi (T k ), H, AljiTi^'WTN) 
27 



(70) 



r 



<k(^ = E c li( r N\\[A mi (^k),H, A Bi (r J )] r ||r JV > 

r 

M (mi)k,(nj)l = ^^(TNWlA^Tk), H, Al^T^fWTN) 
r 

= T,Cl l (T N \\[AUT k ),H } A n3 (T l )] r \\T N ). (71) 
r 

We obtain the matrices N when we replace the double commutators by ordinary commu- 
tators and leave away the Hamiltonian H. As a consequence we have 

N (2) = N (3) = Q (72) 

It is understood that we use the nucleon ground state as defined in Eq. (|40| ) in these 
expressions. Note, however, that with a suitably generalized definition of particle and 
hole states the description remains valid for any ground state with the correct tensor 
structure. 

So far we have suppressed the color degree of freedom. According to the confinement 
property of QCD we demand that the excited states and the ground state are color 
singlets. Then the excitation operator is restricted to color singlet ph pairs. We can 
easily incorporate the color degree of freedom in Eq. (|7T|) by making the replacement 

OjQm * ^ y Q<irP"mc- C^) 

Vo c 

A similar expression holds for the adjoint ph operators. We will display the color quantum 
numbers only when they are explicitly needed. 

In case of a scalar ground state with respect to the group 577(2) j g) SU{2) T the tensor- 
RPA equation (|69| ) can be shown to reduce to an ordinary RPA equation fl30| . 
Given a real interaction, we find that is equal to and that is symmetric in 
case of a time-reversal invariant (in terms of SU(2)j) system, i. e. a scalar ground state. 
For a nonscalar ground state these relations are found not to be fulfilled anymore. How- 
ever, restricting the summation in Eq. ([71]) to T — (0, 0) we recover symmetric expressions 
for the RPA matrices. The symmetry is violated because we need to superimpose terms 

28 



with different T . 

The implication of time-reversal invariance and the resulting structure of the RPA matrix 
is a symmetry in the spectrum of eigenvalues. We find that in this case the eigenvalues 
occur in pairs such that for every uja there is an adjoint — u; A . This way we are guaranteed 
to obtain as many positive as negative eigenvalues, the negative ones are then discarded 
as unphysical. 

On a nonscalar ground state this symmetry is broken. In part II of this paper |34| we will 
find that for the cases under investigation the spectrum is still approximately symmetric 
as long as the eigenvalues are real. 

B. The Basis States 

A genuine feature of the tensor-RPA approach is the overcompleteness of the space of 
basis states defined by Eq. fl67l). A basis state can be labelled by its ph quantum numbers 
(mi, Tk), where m and i denote the quantum numbers of the shells in which particle and 
hole are created and defines the J/T quantum numbers the ph pair is coupled to. A 
novel feature of the tensor RPA is the fact that the ph space for fixed quantum numbers 
Ta of the excited state contains states which are degenerate in energy but carry different 
quantum numbers Tfc. We will distinguish four different classes of basis states: 

A. States with a hole in the Dirac sea and a particle in an empty shell. 

B. States with a hole in the Dirac sea and a particle in the valence shell. 

C. States with a hole in the valence shell and a particle in an empty shell. 

D. States with a hole and a particle in the valence shell. 
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For the moment we do not consider the states of class D, since they require a modified 
RPA scheme. We will come back to this issue later. 

The metric tensor N contains the mutual scalar products of the basis states, so in principle 
we can answer all questions concerning linear dependence and orthogonality by computing 
the relevant matrix elements of N. 

One can easily show that the states of class A are mutually orthogonal and normalized, 
just like in the usual RPA. The states of class B and C do not share that feature, but we 
can derive some general statements about the number of linearly independent states in 
each case. 

We define the degeneracy d of a ph pair as the number of allowed J/T values it can 
be coupled to in order to excite the ground state to a state with quantum numbers 
Ta = (J a, Ta). 

If we restrict ourselves to class B(C) we will always have a particle (hole) in the valence 
shell, the corresponding hole(particle), denoted by r, being in the Dirac sea (in an empty 
shell). In case of the nucleon we obtain for the states of class B and C 



d = < 



1 : j r = Ja ± 1, t r = Ta ± 1 

, jr — J A, t r = Ta ± 1 

2 : { (74) 

jr — J A ± 1) U = T A 



4 : j r — Ja, t r — T A . 

The number of linearly independent basis states for given quantum numbers (mi), which 
we denote by d r , will be less than or equal to d. In general the basis states are of the form 

\(T k ) mi ) = ((alx ai )^x\T N )) rA . (75) 

The ansatz for the RPA operator also contains time reversed ph operators, but the metric 
tensor N can be written in terms of mutual scalar products of the states in Eq. ([75|) only. 
We rewrite this expression by making use of the Racah-recoupling theorem, but we have 
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to make a distinction between the states B and C. 
For the B-states (i. e. states of class B) we obtain 

((a*, x ajfx x 1^))^ 

= ( _ 1) r m +r l -r fc+ i ^ W (^ m r i r A r N , T k T) (a, x (at x \T N }) r f A . (76) 
r 

Effectively this is a unitary transformation from states |(r&) m i) with a coupled ph pair to 
states |r 4g ) with 4 coupled valence quarks. 

In a similar way the C-states are reexpressed in terms of states with two coupled valence 
quarks. Here the situation is slightly more complicated, because we need an explicit 
representation of the nucleon wavefunction. It is not difficult to verify that 

+4„,r,,(«r t „ x a?jM) |0>. (77) 

Only the J/T and color quantum numbers have been displayed. When we insert this 
expression into Eq. (|75|) we get a contribution from each term of the sum. After recoupling, 
the contribution from, e. g. the third term, reads 

I (r k ) mi , 3) = E f f fc H/(r m r i r A r rv , r fc r) 

r 

• [at x (al x (a? w x (a^ x aPJ^f^ |0). (78) 

Corresponding relations are found for the other terms. Bearing in mind that i is a state 
in the valence shell, we can simplify this expression by commuting a? to the right. We 
obtain 



|(r fc ) mi ,3> = ^ff^r^rAr^Tfcr) 

r 
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The first and second term in Eq. ([77]) give similar expressions. This shows that we can 
express the C-states as linear combinations of states |r 2g ) with two coupled valence quarks 
and a third quark in a higher shell. 

As mentioned before, the number of linearly independent B- and C-states is smaller than 
the degeneracy d. The reason for this is the fact that the states |r 4(J ) and |r 2(? ) do not 
all comply with the Pauli principle for any given T^, T2 q - By applying Young-tableaux 
techniques it is found that the value (J=0,T=0) is not allowed for T^ q . For the values 
(0, 1) and (1, 0) are excluded. The derivation of these results is presented in Appendix B. 
For given ph quantum numbers (mi) the forbidden states do not in any case contribute to 
the right hand side of Eq. (|76|) or ([Hj) . Due to the triangle selection rules some of the 6j- 
symbols vanish for certain values T = T'. For this reason the states with the corresponding 
quantum numbers do not appear in the sum over T. If, however, a 6j'-symbol multiplying 
a forbidden state is nonzero the number of linearly independent states d r is lowered by 
one, since the forbidden state has zero norm. An analysis of the triangle selection rules 
for the 6j-symbols reveals a simple relation between d r and the degeneracy d. For the 
B-states we obtain 



whereas for the C-states 



d = 1 
d = 2 
d = 4 

d=l 
d = 2 
d = 4 



10) 



In a practical calculation one would have to reduce the dimension of the naive ph basis 
to d r in each degenerate subspace. The reduced basis need not to be orthogonal, because 
nonorthogonal states are taken care of by the metric tensor N, which contains the mutual 
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scalar products. Alternatively, a Schmidt-orthogonalization procedure might be used. 
Bearing in mind this reduction scheme the basis states of Eq. (|75|) completely span the 
space of physically allowed ph excitations on the nucleon. They are color singlet states 
which belong to SU(2)j <g> SU{2)t multiplets and allow us to treat excitations of the 
valence quarks and excitations of the Dirac sea on the same footing. 
In the nonrelativistic limit the C-states span the space of spin-flavor wavefunctions one 
encounters in the nonrelativistic quark model. There one discards certain spin-flavor 
representations from the very beginning, because they constitute spurious center of mass 
excitations (e. g. the symmetric spin-flavor representation for the lowest f~ , |J state). In 
contrast to that we include the spurious states in our basis and rely on the fact that they 
decouple from the excitation spectrum in a consistent RPA calculation. In a relativis- 
ts framework this is the method of choice, especially when there are broken symmetries 
other than the translational symmetry, like in the case of mesonic modes related to chiral 
symmetry. 



C. The Interaction Matrix Elements 

So far we have only discussed issues which do not depend on the structure of the 
Hamiltonian H . We now focus our attention on the interaction matrix elements and give 
two different prescriptions to compute these. 

We have to compute reduced matrix elements of the double commutators appearing in 
Eq. (|7lD . To do so, we insert the definition of the reduced matrix elements according to 
the Wigner-Eckart theorem. In case of, e. g. MW we get 

(^11^(^)^,4.(^)^1^) = f„ £*^ (82) 

(1 N Y 1 2 U|1 N L z ) 
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This expression is obviously independent of the projection quantum numbers T z , but we 
have to make a choice for T z in order to perform the calculation. In the following we will 
not always display the dependence of all quantities on the projection quantum numbers 
and assume a fixed value of T z in case of doubt. We write the interaction in the general 
form 

H int = \ Y, (PQ\v\RS)aW Q a s a R . (83) 

Z P,Q,R,S 

The sum is taken over the self-consistent single particle states. We use the convention that 
upper case letters denote labels that include the color quantum number (e. g. P — (p, c p )), 
whereas lower case letters label all quantum numbers besides color. 
Usually one splits the interaction into a mean field and a residual part by making use of 
Wick's theorem. The residual part is defined via a normal ordered product with respect 
to the ground state 

N(a ( P a 1 Q a s a R ) ( r N ,T z )- (84) 



With the theorems for normal ordered products |35j it is then easy to compute the RPA 
interaction matrix elements. 

However, the application of Wick's theorem is only possible when we can define a set of 
quasiparticle operators that annihilate the ground state. This is no longer true for the 



nucleon state Eq. ([40]). Even the most general ground state we can write down in the 



framework of our model 

\N')=YCrjT Nl T z ) (85) 

contains correlations beyond the level of Bogolyubov quasiparticles [[H] for any choice of 
the coefficients Cr z . 

To compute the interaction matrix elements one can still assume normal ordering, thus 
defining an approximation that neglects the explicit ground-state correlations which are 
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due to the coupling of the valence quarks. This approximation is solely based on algebraic 
relations for the normal ordering operator and allows us to express the RPA interaction 
matrix elements by the one-body density of the ground state. 
By Eq. (|83|) the residual interaction part of the reduced matrix element becomes 

D R (mi,nj) = l E (PQ\v\RS) ^ - 

•-(r^r^ a\ c a mc ,N(a i P a t Q a s a R )(T NtT ^,al lc ra jc/ \T N ,T Z ). (86) 

In this expression all sums needed for the tensor coupling have been suppressed and only 
the term entering M^> has been displayed. We obtain the analogous expressions for 
M ( 2 ), and if we make suitable interchanges of the quantum numbers m, i, n and 
j. For the evaluation of the normal ordered product we use the theorem |35 

N(d 1 ...d n ^ 1 d n )d n+1 = 

N( o 1 ...d n o n+1 ) + ... + n(6 1 . . .dnd n+1 ) + n{6 1 . . . d n d n+1 ), 

where OjO j = (OiOj) denotes a contraction with respect to jlV, r z ). 

The Oi can be any of the creation or annihilation operators appearing in Eq. (|86|). By 
successive use of the above relation we can express the double commutator in Eq. (]5H) by 
contractions of two operators and arrive at 
1 f 

D R (mi,nj) = - " {(mc)(jc')\v\(ic)(nc')) ■ [p(i) - p(m)][p(j) - p(n)], (88) 

<J \1 AT J- J- 2 U|1 jvl z ) 
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where p is the density matrix as defined in Eq. fl47|) . In this approximation the modi- 
fications with respect to the usual RPA expression are basically the weight factors that 
multiply the ph matrix element. These weight factors account for the fact that the valence 
shell is not fully occupied. 

To compute the mean-field part of the RPA matrix elements, we assume that the one- 
body part of the Hamiltonian, which is the sum of the relativistic kinetic energy and the 
mean-field potential Eq. (pU|), has been diagonalized 
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Hmf = J2 e K a K a K- (89) 

K 

At this stage an average over the projection quantum numbers of the nucleon is assumed to 
identify the mean-field Hamiltonian with the expression discussed in the previous chapter. 
The mean-field contribution to the reduced matrix elements is then defined as 
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DMF{mi 1 ni)=^e K —————--{T Nl T z \ a\ c a mci a\a Kl a\ cl a jd |rjv,r 2 ), (90) 

K \1 N>- 1 «U|1 jv± z) 3 L J 

which can be simplified to give 

D M F(mi, nj) = —————— 6 mn 8 i:j (e m - e*) ■ [p{i) - p{m)\. (91) 

{L N L 1 Z 0\L N L z ) 

Again we find an expression well known from the usual formulation of the RPA multiplied 
by weight factors. 

In order to correctly incorporate the effects that are related to details of the ground-state 
structure, we have to go beyond the normal ordering approximation. Among the basis 
states with lowest energy there will be many B- and C-states, and we expect the lowest 
lying excitations to have a sizeable admixture of these. The B- and C-states are those 
which are sensitive to the structure of the valence shell, as can be seen from the occupation 
factors in Eq. (|88|) . Therefore it is worthwhile to compute the exact interaction matrix 
elements and compare them to Eq. (jB8|). 

For this purpose we do not a priori split the interaction in two parts, but simply compute 
the reduced matrix element of the full interaction 

Dint {mi,nj) = \ £ (PQ\v\RS) ^ 

1 PQRS \L NL L Z \J\L nL z ) 

■g(rAT,r z | a^amcapOQasa^al^ajc/ |rV, T^). (92) 
From this expression we subtract the contributions from the mean-field potential Eq. (|5D|) 



and define the remaining terms to be the residual interaction matrix elements Dr. By 
definition, the mean-field part of the RPA matrix elements is then the same as in the 
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normal ordering approximation. 

For the evaluation of the expectation value in Eq. (|92|) we express the ground state by 
Eq. (|40| ) and commute the annihilation operators in (r7v,r z | to the right until they act 
on the vacuum state |0). We find that there are three different kinds of contributions to 
D R 

D R = D$> + D® + D§>. (93) 

The first one consists of all terms of Eq. ( |5§|) except those weighted with a product of two 
occupation factors g 

D { x\mi,nj) = ^ (r^rroi r^r j ^ mc ^ j ' c '^ I ^ ^ ) 
■ [g(j)0(-i) + g{i)0{-j) - g(m)9(- 3 ) - 9 {n)e(-i) 

—g(j)8(—m) — g(i)8(—n) + g(m)8(—n) + g(n)8(—m) 

+ 18(8(-j)8(-i) - 8(-j)8(-m) - 8(-i)8(-n) + d(-n)8(-m))} . (94) 

The second term only gives a contribution for two ph states with either the same particle 
or the same hole quantum numbers 

1 f 

D^(mi,nj) = — {TnTTz N 01TnTz) GKO " 9(m) + g(j) - g{n)) 

• [Smn ((jc)(kc')\v\(ic)(kc')) - ((mc) (kc') \v\ (nc) (kc'))} . (95) 

k 

The third term cannot be expressed in terms of the occupation factors g, but contains 
the coupling coefficients T explicitly 



D R (mi,nj) (rjvrrz0|r ^ } 

P,Q,R,S 

~\~Sij&qm,G r s,np ~\- 5 mn S r iGpqj s 4:S sn Sq m Gipj r ^$pj $riG n q t rns\ • (96) 

The projection on color- singlet configurations is ensured by 
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C pq ,rs — 8 CpC 5 CqC >(5 CrC >5 CaC — 5 CrC 5 CgC '). (97) 

c,c' = l 

c=£c> 

The color operator C pq ^ rs represents those contributions from the two-body interaction 
a^pC^QasciR where a creation and an annihilation operator are coupled to an intermediate 
color-octet. Note the + sign in front of the exchange term in Eq. fl9"E[), which is due to 
the antisymmetry of C pq ^ s under exchange of r and s. 

The structure of the valence shell is reflected in a matrix G, which is defined as 

\ Er z Tr zp r zq r z Tr zr r zs r z ■ P,q,r,s e valence shell 

(jpq,rs — < . (98) 

I : else 

(2) (3) 

Through terms D R and D R new types of matrix elements are introduced to the RPA 
theory. Eq. (|95|) describes the scattering of particles and holes by the mean field of the 
valence quarks. In addition to scattering processes (i. e. processes in which particles or 
holes keep their identity) Eq. (|96| ) also describes the transition of particles and holes into 
valence states. Similar structures are found in QRPA, where ph and pp/hh matrix ele- 
ments contribute. 



D. Description of (| , |) states 

With the formalism presented so far it is not possible to study excitations of the 
nucleon in the (| + , |) channel. The three valence quarks can as well be coupled to the 
quantum numbers of the A(1232), which can be achieved by acting with a ph operator of 
special structure on the nucleon ground state 

|r A ) = (aUt) x \r N )) r \ (99) 

where T = (1, 1). In this expression both the particle m and the hole i are in the valence 
shell. If the mean-field level of the theory is defined in the sense of the previous section 
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(with an averaged ground-state density) such a state is clearly energetically degenerate 
with the nucleon. For the ph operator in Eq. (|99|) we have 

aUt) = A mi (r). (ioo) 

The adjoint of an RPA excitation operator including contributions from the recoupling 
operator Eq. ( |9"9"D will in general not fulfill the relation 

O xTk \T N )^0, (101) 



since the ph annihilation operators taken from Eq. (|99l) are of the same structure as the 
creation operators. Hence, there is no reason to believe that the Y- Amplitudes multiplying 
the recoupling contribution are small. However, in the derivation of the RPA equation 
we have assumed that Eq. |101| is fulfilled. 

From investigations of hyperfine splitting in the MIT-bag model it is strongly suggested 
that the wavefunction of the A(1232) has a dominant contribution from the recoupling 



configuration ||36|| . Thus an RPA approach in a ph space without the state Eq. fl9H|) can 
at best account for the higher lying resonances in the (| + , |) channel. 
In the following derivation we neglect the tensor coupling of operators and states in order 
to point out more clearly the special features of the equations for the A-channel. A 
straightforeward way to account for admixtures of both types of configurations is to write 
down the full Hamiltonian in a model space containing the recoupled ground state as well 
as the nondegenerate ph excitations. For this purpose we write an excited state as 

|xr A ) = |A„> + |A fc ), (102) 

where x denote additional quantum numbers. The projection on the recoupling configu- 
ration is given by 

|A„) = Z v \v), (103) 
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where \v) is the properly normalized state of Eq. (0). The projection on the space of ph 
states (of nonzero energy) is given by 

\A k )=Y f Z ki \k i ), (104) 

i 

where 

\h)=r,l\Y N ) (105) 

and r)l. are ph operators including the states of class A,B and C. Again we assume proper 
normalization of the basis states. All quantum numbers, including angular momentum 
and isospin, have been absorbed in the label fcj. Clearly, the two subspaces are orthogonal, 
i. e. (ki\v) = Vz. Configuration mixing is contained in the amplitudes Z v , Z ki , where 

\Z V \ 2 + Y.\Zk z \ 2 = 1. (106) 

i 

To proceed, we split off a residual interaction V from the Hamiltonian, in the same way 
as this was done in the previous section. We assume that the mean-field part of the 
Hamiltonian is diagonal in the basis we have chosen. The projection of the eigenvalue 
equation on the different subspaces is given by 

(H vv — E)Z V + VvkiZk, = 

i 

- 5 k^E)Z k . + V kiV Z v = 0. (107) 

3 

Clearly H^. and V v ^ are matrices according to the dimensionality of the ph space. The 
matrix elements are in general divergent due to contributions from the Dirac sea. A 
prescription to regularize these matrix elements in agreement with the cut-off scheme for 



the mean-field Hamiltonian has been given in |p7 |. When we solve the first equation for 



Z v and insert the result into the second, we obtain 

J2( H k ikj + *W£) - S kikj E)Z k . = 0, (10* 

3 
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where a self-energy 





(109) 



accounts for the coupling to the valence configuration. G v is the projection of the many- 
body propagator and is given by 



Thus in the ph channel we obtain an equation where the residual interaction is supple- 
mented by a self-energy term, compared to the conventional formulation without coupling 
to the valence configuration. We immediately find that the self energy will give a con- 
siderable contribution to Eq. (|108|) for energies close to the pole value H aa . Since Eq. 
( |107|) is a diagonalization problem for a hermitian operator, the lowest eigenvalue will be 
located below H aa . 

Denoting the energy-expectation value of the nucleon ground state Eq. (|40| ) by Eq, we 
obtain for the N-A mass splitting 



This quantity must be compared to the phenomenological value for the mass splitting of 
about 300MeV. In the last section it was argued that H aa — E should already be a rather 
small quantity for reasons of consistency with the assumptions entering the derivation of 
the RPA equation. Thus it might be difficult to reproduce the observed value for the 
mass splitting on the level of Eq. (|111|) . 

A solution to this problem could be provided when we treat the nucleon ground state on 
the same basis as the recoupling configuration. Due to the averaging of the ground-state 
density over projection quantum numbers in the course which information about the de- 
tailed structure of the nucleon is lost we can lower the nucleon energy Eq by infinitesimal 
lp — lh admixtures. In the usual Hartree-Fock description such admixtures are excluded 



AE N/ \ < H aa — 



Eq- 



(111) 
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by the self-consistency condition. The derivation presented in this section effectively takes 
into account lp — lh admixtures in the (| + , |) channel to the recoupling configuration 
| a). In a similar way we have to account for lp — lh admixtures to the nucleon in the 
|) channel. This is required by reasons of consistency when comparing the (^ + , ^) 
and the (| + , §) channel. By such admixtures the energy of the nucleon will be lowered, 
thus the mass splitting is increased with respect to the value obtained from configuration 
mixing in the A-channel. 

In addition to that we have to take into account another important point. So far the spec- 
trum of (| + , |) states was derived from a diagonalization problem of the full Hamiltonian 
in a model space. It is important to note that the usual RPA scheme is not equivalent 
to such a diagonalization problem. However, even when we choose excitation operators 
ry k . of RPA form in Eq. (^) we still obtain a diagonalization problem, since the contribu- 
tion of the time reversed operators vanishes. In the subspace of ph states we recover the 
well known RPA structure when we replace the operator products in Eq. (|107|) by double 
commutators 

(r N \r) ki Hril.\r N ) = (r^i^.ir^Ijir^ + ^r^^^iir^). (112) 

Note that within the space of ph operators we have 

v ki \r N ) = o, (113) 

which is an approximate identity when RPA operators are acting on a HF ground state. 



This justifies the replacement in Eq. (|112|) . Note that we are not allowed to introduce 
double commutators in the full space of excitation operators, since the recoupling operator 
does not annihilate the ground state. Abbreviating the recoupling operator of Eq. (|99|) 
with 7)1, we are, however, allowed to subtract the ground state energy 

(r N \ Vv H V l\r N ) = (r N \ Vv [H, v l}\r N ) + E (r N \ VvV l\r N ). (iu) 
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Together with Eq. ( |112| ) we thus obtain an eigenvalue equation for the excitation energy 
c<ja = Ea — Eq. It is easy to show that the replacements of Eq. (|112|) and ( |114| ) preserve 
the symmetry of the eigenvalue problem (i. e. hermiticity) if Eq. ( |113| ) is fulfilled exactly, 
which is the case on Tamm-Dancoff level. For excitation operators of RPA type we obtain, 
as usual, a diagonalization problem for a nonhermitian matrix. Wheras in case of a scalar 
system such an RPA problem can in general be reformulated as an eigenvalue equation for 
a hermitian matrix pl| , the tensor RPA scheme does not allow for such a reformulation. 
This fact is connected to the asymmetry of the spectrum under inversion of the energy 
axis. 

In such a more general formulation Eq. ( |111|) looses its validity. The upper bound for 
the N-A mass splitting is based on the Rayleigh-Ritz variational principle for hermitian 
matrices. Thus for the description of (| + , |) states it seems to be crucial to go beyond 
the Tamm-Dancoff level and introduce commutators instead of operator products. 
On RPA level there is no a priori estimate for the mass of the A(1232) which stands in 
contradiction to the consistency requirements set by the Hartree-Fock theory. Only in a 
numerical calculation for a specific model it can be decided if such a RPA scheme allows 
for a realistic description of the (| + , |) spectrum. 



VI. SUMMARY, DISCUSSION AND CONCLUSION 

The theory presented in this paper is intended to provide a consistent approach to 
many-body effects in the ground state and excited states of the nucleon. By using a 
fieldtheoretical formulation essential features of low-energy hadron physics are accounted 
for. The theory is based on an effective two-body quark-quark interaction which must be 
determined empirically. The approach describes baryon spectroscopy by an expansion in 
terms of many-body correlation functions, where at the present stage RPA contributions 
were investigated. 
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The nucleon is considered as a system of three valence quarks which interact with the 
vacuum. The ground state of the full system, including valence and vacuum parts, was de- 
scribed with Hartree-Fock methods which clearly represents a nonperturbative approach. 
By introducing a symmetry conserving regularization of the mean-field potential a self- 
consistent RPA scheme including excitations from the Dirac sea was defined. 
A central issue of the model is that states of good angular momentum J and isospin 
T are used. This we consider as an important advantage over other approaches which 
are based on a grand spin description, as e. g. the hedgehog soliton |7|,§]. However, it 
was discussed that this leads to mean-field self-energies violating rotational and isospin 
invariance. Although many-body theory provides in principle methods to treat such a 
problem, e. g. in nuclear structure physics ||21|| , we have chosen a somewhat simplified 
approach by averaging over projection quantum numbers. This led to a one-body density 
matrix which is the same for the (| + , |) nucleon and the (| + , |) A-like configuration. A 
mean-field Hamiltonian defined with the averaged density allows to extract single particle 
states with good J/T. From these a nucleon state can be built which contains important 
details of the many-body structure related to angular momentum and isospin coupling. 
Excited states of the nucleon were described by RPA-methods. Since we are dealing with 
a nonscalar ground state the theory must account for the mixing of ph configurations with 
different J/T quantum numbers in an excited state. Appropriate methods are provided 
by the tensor RPA which originally was introduced in nuclear structure physics by Rowe 
et al. |13|. Although technically rather involved, the tensor RPA method relies on a sim- 
ilar representation of excited states as collective lp — lh configurations on a mean-field 
ground state as the standard RPA scheme. In the present case, this includes transitions 
from the valence sector and the Dirac-sea into positive energy states. The application 
of RPA theory to the vacuum sector is well established for the NJL-model 0. As an 
extension of such techniques we have made a first attempt to describe a baryon with good 



44 



J/T quantum numbers in the ground state and in the excited states. 
The description of excited states provides additional information on the residual interac- 
tion beyond those properties which contribute to the ground state. It was pointed out 
that the regularization of the Hamiltonian led to a modification of the residual interaction. 
This reflects the fact that a projection onto a model space is used which includes vacuum 
states only up to a certain cut-off energy. For practical purposes it might be preferable 
to use a parametrization of the residual interaction and account for the modifications 
implied by the cut-off scheme by a readjustment of the parameters. 
It was also pointed out that the RPA theory allows to obtain information on the ground- 
state structure beyond mean-field level. On the one hand there will be lp — lh admixtures 
which are due to the averaging prescription for the ground-state density. On the other 
hand an RPA ansatz includes implicit 2p — 2h correlations. An improved ground-state 
wavefunction including such contributions will allow to recover the coupling to pion de- 
grees of freedom which are excluded on mean-field level for reasons of symmetry. Further- 
more ground-state correlations will lead to a depolarization of the valence-quark core and 
allow to study the spin transfer to, e. g. excitations of the Dirac sea. Such admixtures 
have already been investigated in [52] using a somewhat simplified RPA scheme. 



APPENDIX A 

In this appendix we derive an operator equation for the effective ph interaction which 
has to be used in the Dirac-RPA equation. 

By definition, the mean-field Hamiltonian depends on p only via the the regularized 
density p' 

9 -H' MF = J-H MF ( P ') = ( ^-H MF ) Ml. (115) 



Opkk' Opkk* \OPw J op. 



kk' 
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We have adopted the convention to sum over repeated indices. The regularized density 
depends on p also via the regularization function R 

dp'iv _dR K dp* , Rn dIU>v , . 

OPkk' OPkk' OPkk' OPkk' 

We arrive at 



~H'mf — [ a i Hmf) RlkRlk' 



dpkk' \ d P'u> 

+ Tri-J^MF Pi; 



r\ / IV1 r I f Vt c\ 

PPlV / OPkk' 



dRii dRi'i' 
-tivv + tin- 



d(H' MF ) mm , u d(H' M 



(117) 



' mm' mm' J 

The effective interaction dH' MF /dp appears on both sides of the equation. Expression 
( |117|) has to be evaluated at po, the stationary solution of the TDHF equation. For the 
remaining parts of this section we adopt the convention that single particle labels refer to 
the basis in which po is diagonal. 

Denoting the eigenvalues of H' MF by e^, we find that the first term on the right hand side 
of Eq. (|117| ) is given by 

( d 



-H 



MF 



R[e k ]R[e k ,} (118) 



where all density-dependent quantities have been evaluated at p . In the self-consistent 
basis the density matrix acquires the form 

(po) kk , = 5 k k>0(k F -k), (119) 

for states normalized to unity. The symbolic ^-function indicates that k has to be below 
the Fermi surface which is defined by hp. In an infinite system hp denotes the Fermi 
momentum. 

To evaluate the second term in Eq. ( p. 171 ), we write the regularization function as 



R[x] = a n x n . (120) 

n=0 

For an arbitrary operator A we have 
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OR 



dA Ul 



oo 

--H^ny-^— 

mm' n=0 \ u - rL rnm 



PlA 



A ■ ■ A ■ 



z *i"fTl • • • Ai n _ii + . . . + Ai^Ai^ ... 



u - rl -mm' u - ri -mm' / 

where the Einstein convention for repeated indices is used. Evaluating this expression for 



A = h 'mf g ives 



OR, 



X X \ 1 „ / n— 1 i ,n— 2 , , ,n— 2, , , 

— OlmOim' a n\^m' + e m e m' + • • • + € m e m' + e 

PO «=0 



71— 1\ 

m / 



:i22i 



The expression in brackets can be cast into the form of a finite geometrical series. Then 
the derivative of the cut-off function simplifies to 



dR, 



&lm fiim' 



po €m ' tm \n=0 



n=0 / 



(123) 



d\H'j^p) mm i 

On the right hand side we recover the defining power series for the cut-off function R and 
obtain 



dR, 



dm 



M F ) mm' 



^Im^im' 



R{e m >] — R[e r , 



(124) 



Po 



With that we finally obtain an operator equation from which the effective interaction can 
be determined 



d 



-m 



dpkk' 
dH MF 



MF 



Pa 



f IF R[e k }R[€k> 



(125) 



+ 



dp'w 



d { h 'mf) w 



Pa 



dp 



kk' 



Pa 



(^-o^hl + ^^-O^M)- 



APPENDIX B 

The structure of 3-quark wavefunctions has been discussed in detail in the literature, 
e. g. in the book by Close f3"8|| . In this appendix we investigate the group representations 
which can build from the 2- and 4-quark wavefunctions we encounter when we recouple 
the RPA basis states. 
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First we discuss the B-states, which contain a hole in the Dirac sea and four quarks in 
the valence shell. 

A single particle state is characterized by its quantum numbers in Dirac, isospin and color 
space. We assume that the valence shell is a \k\ = 1 state, so that the Dirac-space part of 
the wavefunction is uniquely determined by specifying the angular-momentum projection 
quantum number. We can simply imagine having spin | particles with additional isospin 
and color degrees of freedom, so we will label the Dirac-space representations with a 
subscript J, denoting the total angular momentum. The isospin and color representations 
carry subscripts T and C, respectively. 

We have to construct irreducible representations of the group SU (2) j ® SU (2) T ® SU (3) c ■ 
In order to be allowed by the Pauli principle, these states have to form antisymmetric 
representations of the symmetric group. The space of allowed 4-quark wavefunctions is 
further restricted to those combinations which can be reached by coupling one quark to 
a (J=|, r=|) color singlet 3-quark state. 

According to the rules for outer products of Young tableaux the only relevant 

SU(2) diagrams are (3, 1) and (2,2) with dimensionalities 

I 3 : (3,1) 

dsu{2) = { (126) 
| 1 : (2 2 ) 

When we take the inner product of a J with a T representation we can classify the resulting 
representation by a SU(4)jt Young tableaux. We could now deduce the SU(2)j®SU(2)t 
content of the £77(4) representations by using the Clebsch-Gordan coefficients of the sym- 
metric group ||40|| . In many cases this is not neccessary, since one often can conclude on 
which product representations of lower dimensional unitary groups are contained in a 
Young tableaux by considering the symmetry of the diagrams and matching the dimen- 
sionalities of the representations. 

One finds that the product representations of the relevant J and T diagrams are contained 
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in SU (4) j T diagrams, according to 

SU(2)j x SU(2) T 



SU(4) 



JT 



(3,l)x(3,l) (4),(3,1),(2 2 ),(2,1 2 ) 

(3,l)x(2 2 ) (3,1),(2,1 2 ) (127) 

(2 2 )x(3,l) (3,1),(2,1 2 ) 

(2 2 )x(2 2 ) (4),(2 2 ),(1 4 ) 
The SU(3)c representation is fixed to be (2, l 2 ), since the three valence quarks of the 
nucleon form a color singlet. Coupling this C-representation on the SU(4) diagrams one 
can embed the resulting product representations in multiplets of the group S77(12) jtc, 
according to 



SU{A) JT x SU(3) C 



(4) x (2, l 2 ) 
(3,l)x(2,l 2 ) 
(2 2 ) x (2,1 2 ) 
(2,l 2 )x(2,l 2 ) 
fl 4 ) x (2, l 2 ) 



SU{\2) 



JTC 



(128) 



(2, 1 2 ) 

(1 4 ),(2, 1 2 ),(2 2 ),(3,1) 
(3,1),(2,1 2 ) 
(4),(3,1),(2 2 ),(2,1 2 ) 
(3,1) 

The 577(12) representation is restricted to the completely antisymmetric (l 4 ). Therefore 
the only allowed SU{A) diagram is (3, 1). From the spin-singlet isospin-singlet diagrams 
we cannot construct this SU (4) representation, so a configuration of four valence quarks 
with ( J=0,T=0) is excluded by the Pauli principle. 

For the states of class C we can apply quite similar arguments based on a configuration 
of two valence quarks. To obtain the possible 2-quark Young tableaux, we again label the 
Dirac space diagrams with a subscript J, since all other quantum numbers are the same. 
Then the dimensionalities of the Dirac and isospin diagrams are 

3 : (2) 
1 : (I 2 ) ' 



d 



SU{2) 



(129) 
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These representations can be coupled according to 



SU(2)j x SU(2) T 


SU(A) JT 


(2) x (2) 


(2) 


(2) x (I 2 ) 


(I 2 ) 


(I 2 ) x (2) 


(I 2 ) 


(I 2 ) x (l 2 ) 


(2) 



When we take the outer product of the 2-quark SU (4) diagrams with a third particle in 
a higher shell, we obtain 

(2)0(1) = (3) + (2,1) 

(1 2 )®(1) = (2,1) + (1 3 ). (131) 

In color space the three quarks form a completely antisymmetric representation, so the 
residual wavefunction has to be symmetric under exchange of particles. Only the sym- 
metric SU (4) representation for two quarks can be coupled to a symmetric 3-quark rep- 
resentation, so the (J=0,T=1) and the (J=1,T=0) 2-quark configurations are forbidden 
by the Pauli principle. 
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